CCFM prediction on forward jets and F2'. 
parton level predictions and a new hadron level Monte 

Carlo generator CASCADE 

^ ; H. Jung 

Physics Department, Lund University, Box 118, S-221 00 Lund, Sweden 



^3 



OO 



Abstract: A solution of the CCFM equation for a description of both the structure 



^ ■ function F2 and the cross section of forward jet production as measured by the 



HERA experiments is obtained on the basis of the parton level Monte Carlo program 

\ Smallx. The treatment of the non - Sudakov form factor and the so - called 

"consistency constraint" are discussed. Following this a backward evolution scheme 

^ ■ according to CCFM is developed which then is used to construct an efficient hadron 

■ level Monte Carlo program Cascade. The results from the forward evolution and 

. the backward evolution Monte Carlos are compared and found to be consistent. 

00 ■ 
O 

a^ ■ 1 Introduction 

^ The parton evolution at small values of x is believed to be best described by the CCFM evolution 

Q-i! equation [ |I|, ^, ^ , which for a; ^ is equivalent to the BFKL evolution equation [ §, §, 0] 

r-| , and for large x reproduce the standard DGLAP equations. The CCFM evolution equation 

\ takes coherence effects of the radiated gluons into account via angular ordering. 

^ ' Already in 1992 the CCFM evolution equation was implemented in the parton level Monte 

5-1 Carlo program, Smallx [ 0, using a forward evolution scheme. In 1997 the Linked Dipole 

■ - - ' Chain [ |10|, |ll], |12|, |13| and the corresponding LDC hadron level Monte Carlo program were 



developed based on a reformulation of the original CCFM equation. Predictions of the CCFM 



equation for hadronic final state properties were studied in [ 14], paying special attention to 
non-leading effects. Common to all three approaches was the difficulty to describe the structure 
function F2 and the cross section of forward jet production in deep inelastic scattering at the 
same time [ |rU|, |n|, 0, [T|, |g, |T§ . 



This article is divided into two parts: In the first part I shall describe briefiy the imple- 
mentation of the CCFM evolution equation into the program Smallx [ ||, 3, and discuss the 
treatment of the non-Sudakov form factor A„s as well as the effects of the so-called "consistency 
constraint" , which was found to be necessary to include non-leading contributions to the BFKL 
equation [ . Then I demonstrate that a proper treatment of the non-Sudakov form factor is 
essential to achieve a good description of F2 and the forward jet data at the same time. 
In the second part I describe a model for a backward evolution scheme according to the CCFM 
equation, which forms the basis of the new hadron level Monte Carlo generator Cascade. One 
ingredient in this approach is the determination of the unintegrated gluon density function, 
which is obtained from the forward evolution in Smallx. As a check of consistency I compare 
results obtained from Smallx with those from Cascade and show nice agreement. 
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Figure 1: Kinematic variables for multi-gluon emission. The t-channel gluon four - vectors are 
given by ki and the gluons emitted in the initial state cascade have four - vectors pi. 

2 Forward evolution: CCFM and Smallx 

The implementation of the CCFM [ |I], 0, ^, Q parton evolution in the forward evolution Monte 
Carlo program Smallx is described in detail in [ ||, § . Here I only concentrate on the basic 
ideas and discuss the treatment of the non-Sudakov form factor. 

The initial state gluon radiation is sketched in Fig. |I| together with a definition of some 
variables used in the following. According to the CCFM evolution equation, the emission of 
partons during the initial cascade is only allowed in an angular ordered phase space region. 
The maximum allowed angle S is defined by the hard scattering quark box, which connects 
the gluon to the virtual photon. In terms of Sudakov variables the quark pair momentum is 
written as: 

Pq + Pq = Y{Pp + Spe) + Qt (l) 

where Pp and pe are the proton and electron momenta, respectively and Qt is the transverse 
momentum of the quark pair. Similarly, the momenta Pi of the gluons emitted during the initial 
state cascade are given by (here treated massless): 

Pi = ViiPp + iiPe) +Pti , 6 = (2) 

^Vi 

with Hi = {1 — Zi)xi^i and Xj = and s = {pp +PeY being the total electron proton center 

of mass energy. The variable is connected to the angle of the emitted gluon with respect to 
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the incoming proton and Xi and i/i are the momentum fractions of the exchanged and emitted 
gluons, while Zi is the momentum fraction in the branching (i — 1) i and pu is the transverse 
momentum of the emitted gluon. 

The angular ordered region is then specified by: 

Co < 6 < ■ • • < a < 2 (3) 

which becomes: 

Zi^iqu-i < qti (4) 
when using the rescaled transverse momenta qu of the emitted gluons defined by: 

qti = Xi_i^/s^i = (5) 



In Smallx, the initial state gluon cascade is generated in a forward evolution approach 
from a initial distribution of the kt unintegrated gluon distribution according to: 

xG,{x, kl,) = N-{l-xf- exp (-fc^o/fc^^) 

where is a normalization constant. The exponential factor in eq.@ is a gaussian distribution 
which specifies the distribution of the initial transverse momenta kto- In the following we set 
/cq = 1 GeV^, corresponding to a gaussian width of a = fco/v^ = 0.7 GeV. The input gluon 
distribution needs to be adjusted to fit existing data, but it turns out that the small x behavior 
of the structure function F2 is rather insensitive to the actual choice of xGq{x^ kj) and only the 
normalization acts clS cl free parameter. 

The initial state branching of a gluon ki into another virtual (t-channel) gluon fcj+i and a final 
gluon pi+i (treated on mass shell) is generated successively from the inital gluon distribution at 
a starting scale Qo- The probability for successive branchings to occur is given by the CCFM 
splitting function [ |T], |, |, |: 

dPi = PHz,, ql kl) ■ AJz,^ ■ e(gi - z,g,_i) • 6(1 - - e,) (7) 

Tiqi 

with qi = Pti/{1 — Zi) being the rescaled transverse momentum of the emitted gluon i. The 
fractional energy of the exchanged gluon i is given by Xi and the energy transfer between the 
exchanged gluon i — 1 and i is given by Zi = Xi/xi-i. A collinear cutoff = Qo/qi is introduced 
to avoid the 1/(1 — singularity. The Sudakov form factor is given by: 

A,{qi, Ziq,.i) = exp - / / dz^^ ^ 8 

\ q^ Jo 1- z J 

with a, = = The Sudakov form factor cancels at an inclusive level against the 

1/(1 — 2;) collinear singularity of the splitting function. Coherence effects are taken into account 
by angular ordering qi > Zi^iqi^i given by the first 9 function in eq.(0). The distribution of 
successive values of is controlled by this Sudakov form factor. The cascade continues until 
reaches Pmax-, which is given by the partons from the hard scattering matrix element. 
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The gluon splitting function is given by: 



1 - Zi 



7- n"^ 

nsy^tj Hi } i^ti) 



(9) 



with kti being the transverse momentum of the exchanged gluon i and the non-Sudakov form 
factor Ans being defined as: 



log = -asik'^i) J ^ J %^Q{kti - q)Q{q - z'qt. 



q' 



(10) 



The difference to the DGLAP splitting function for gluons is the appearance of the non-Sudakov 
form factor Ans, which screens the 1/z singularity in eq.(^. The finite terms in the splitting 
function are neglected, since they are not obtained in CCFM at the leading infrared accuracy [ 
^ p. 72]. In the region > Ziqf, Ans can be written as: 



log A„, = -as{kl) log - I log 



(11) 



where the limits of the z' integral in eq. ([lOD were set to Zi < z' < 1. However, the functions 
in eq.(p!0|) give further constraints on the upper limit of the z' integral^: 



kf, 



(12) 



Zi < z' < Zq = max | 1, — 
In the region Zi > ku/qi, the z' integral becomes zero. Putting all this together the non-Sudakov 

as{kl) 



form factor can be expressed as [ |T^ 

log = 
where 



2n 



log — log 



^ti 
^^oZiqf 



(13) 



Zq 




if kti/qi > 1 
if Zi < kti/qi < 1 
if kti/qi < Zi 

which means that in the region kti/qi < Zi we have Ans = 1, giving no suppression at all. 

The use of eq.(^T]) in the region kfi < Ziqf = ZiPti/{l — Zi)^ would result in a change of 
the sign, and therefore instead of giving a suppression form factor would give an enhancement, 
which is in contradiction to a form factor. The constraint]^ kf^ > Ziqf is often referred to as the 
"consistency constraint" [ [16|. 

am very grateful to J. Kwiecinski for the explanation of these limits 
^ The constraint: < Qf/zi was already derived in [ |l8[ eq.(2.11)] and in [ |l^ in a frame where parton 

i is along the z axis. In [ ^ this was used to derive a similar condition on the transverse momenta (see [ ^, 
eq.(7.23)]) fcj > zq^, when the transverse momenta are the dominant contribution to the virtualities. However, 
in a different frame (the lab frame used in CCFM for example), the numerical value of z might be different, 
this constraint might be violated. 
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In the original version of Smallx [ H, 0] the sphtting function was defined as: 



which also is valid over the full phase space, but gives only very weak suppression of the small 
z splittings, and therefore overshoots the measured structure function F2 at small values of x. 

The Monte Carlo program Smallx [ ||, § has been modified to include the non-Sudakov 
form factor according to eq.(|T3|) and kf was used as the scale in also in the matrix element. 
To avoid problems at small fc^, a upper limit on was introduced such that as{k^j) < 0.6. 
This modified version of Smallx is called Smmod in the following. 

A few words are needed concerning the "consistency constraint", which was introduced to 
account for next-to-leading effects in the BFKL equation. Including this constraint it was found 
in [ that about 70% of the full next-to-leading order corrections to the BFKL equation are 
simulated. However, in LO BFKL the true kinematics of the branchings are neglected. They are 
included only in NLO therefore act as next-to-leading order effects. Therefore this constraint 
is often also called "kinematic constraint". In the CCFM equation energy and momentum 
conservation is already included at LO by the full treatment of the radiated gluons, and it is 
not clear, whether the arguments coming from BFKL also apply to CCFM. 

In the following the effects of the "consistency constraint" are studied in more detail. In 
Fig. 1^ the structure function ^2(0;, Q^) as calculated from Smmod (including the non-Sudakov 
form factor according to eq.(|l^)) is shown. The solid curve is without the "consistency con- 
straint", the dashed curve shows the result after applying the "consistency constraint". The 
collinear cutoff Qo was set to Qo = 1.1(0.85) GeV (without (with) "consistency constraint", 
respectively). The values of Qq have been chosen such, that a reasonable description of F2 
is obtained. In Fig. ^ the cross section of forward jet production and the measurements of [ 
2^, ^ ^ are shown as a function of x (Fig. ^ja.-c.) and as a function of E^/Q"^ (Fig. 0d.), using 
the same parameter setting as for the calculation of F2. This plots have been produced with the 
HZTOOL analysis package [ 2l|. Nice agreement with the data is observed for the case without 



the "consistency constraint" . When the "consistency constraint" is applied, the prediction falls 
below the data. One has to note that in all cases the scale for is fc^ . Comparing this to the 



results in [ the same trend is found: applying the "consistency constraint" and using kf as 



the scale in as, the prediction of the cross section for forward jet production falls below the 
data, whereas there is essentially no effect seen in the prediction of the structure function F2. 
It is interesting to note that the BFKL result of [ (dashed curve in Fig. 4 of [ |TB|) seems to 
be identical to the dashed line in Fig. ^ obtained from Smmod. 



3 Backward evolution: CCFM and CASCADE 

The forward evolution procedure as implemented in Smallx is a direct way of solving the 
CCFM evolution equation including the correct treatment of the kinematics in each branching. 
However the forward evolution is rather time consuming, since in each branching a weight factor 
is associated, and only after the initial state cascade has been generated completely, it can be 
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Figure 2: The structure function F2{x, Q"^) compared to HI data IW^- The solid (dashed) line is 
the prediction of the Smmod Monte Carlo without (with) applying the "consistency constraint" 
(c.c). 



decided whether the kinematics allow to generate the hard scattering process. Quite often, a 
complete event has to be rejected. 

A more efficient procedure to be used in a full hadron level Monte Carlo is a backward 



evolution scheme, which is used in standard Monte Carlo programs [ ^ using a DGLAP 
type parton cascade. The idea is to ffist generate the hard scattering process with the initial 
parton momenta distributed according to the parton distribution functions. This involves in 
general only a fixed number of degrees of freedom, and the hard scattering process can be 
generated quite efficiently. Then, the initial state cascade is generated by going backwards 
from the hard scattering process towards the beam particles. In a DGLAP type cascade the 
evolution (ordering) is done usually in the virtualities of the exchanged t-channel partons. 

According to the CCFM equation the probability of finding a gluon in the proton depends 
on three variables, the momentum fraction x, the transverse momentum squared k'j of the 
exchanged gluons and the maximum angle allowed for any emission Pmax = a;„_i\/sH. This 
probability distribution as a function of x, and Pmax is obtained from a solution of the 
CCFM equation, here taken from the Monte Carlo procedure used in Smallx and Smmod as 
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Figure 3: a. — c. The cross section for forward jet production as a function of x, for different 
cuts in pt compared to HI data f\2l\] (a. — b.) and compared to ZEUS data /[^/ (c). d. The cross 
section for forward jet production as a function ofE^/Q^ compared to /[^/. The solid (dashed) 
line is the prediction of the Smmod Monte Carlo without (with) applying the "consistency 
constraint" (c.c). 



described in the previous chapter. For practical reasons the gluon density is obtained in a 
50 X 50 X 20 grid in logx, log kt and log Pmax, and then a hnear interpolation method is used 
to obtain the gluon density at values in between the grid points. 
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However there is a problem in the normahzation of the gluon density: In Smallx and 
Smmod the normahzation is done via the momentum sum rule for all gluons below a given 
kt for all Pmax- But this Pmax depends on the kinematics of the process involved. In order 
to obtain a gluon distribution, which is universal, the momentum sum rule is applied for each 
Pmax separately. Therefore differences between Smmod and the cross section calculated in the 
following are not surprising. 

To fully simulate the hadronic final state of a small x process including the CCFM parton 
evolution, a new hadron level Monte Carlo program. Cascade, has been developed. It consists 
of three steps: 

• First the hard scattering process is generated: 

cr = J dkldXgA{xg, k1, Pmax)(^il*9* qq) (15) 

where the off-shell matrix elements (given in [ ^ p. 178 ff]) are used, with the gluon 



momentum (in Sudakov representation): 

k = XgPp + XgPe + kt XgPp + kf (16) 

The gluon virtuality k^ is then k"^ kf. 

The initial state cascade is generated according to CCFM in a backward evolution ap- 
proach (described in the next section). 

The hadronization is performed using the Lund string fragmentation implemented in 
JETSET |28|, H. 



At the present state, final state gluon radiation from both the quarks and the gluon pi are 
neglected. 

However in the backward evolution there is one difficulty: The gluon virtuality enters in the 
hard scattering process and also influences the kinematics of the produced quarks and therefore 
the maximum angle allowed for any further emission in the initial state cascade. This virtuality 
is only known after the whole cascade has been generated, since it depends on the history of 
the gluon evolution. In the evolution equations itself it does not enter, since there only the 
longitudinal energy fractions Zi and the transverse momenta are involved. This problem can 
only be approximatively overcome by using k'^ = k^/{l — Xg) for the virtuality which is correct 
in the case of no further gluon emission in the initial state. 

The Monte Carlo program Cascade can be used to generate unweighted full hadron level 
events, including initial state parton evolution according to the CCFM equation and the off - 
shell matrix elements for the hard scattering process. It is applicable both for photo production 
of heavy quarks as well as for deep inelastic scattering. The typical time for generating one event 
is ~ 0.03 sec, which is similar to the time needed by standard Monte Carlo event generators such 
as LEPTO [ ^ or PYTHIA [ In standard Monte Carlo generators [ |30|, uncertainties 



concerning the scale to be used in the the structure functions and in a^, and also the upper 
virtuality used in the parton shower approach are typical. As in Smmod and Smallx, in 
Cascade, there are essentially no such free parameters left since everything is fixed by the 
requirement to describe the inclusive structure function F2. 
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3.1 Backward evolution formalism 



The CCFM evolution equation can be written [ ^, [1^ in a integral form as: 

— / — ^Q{Pmax - zq)As{Pmax,zq)P{z,q,kt)A{-,k[,q) 

(17) 

with k[ = \kt + {1 — z)q\ and with Pmax being the upper scale for the last angle of the emission: 
Pmax > Znqn, In > -z^n-iQ'n-i, li > -^o^o = Qo- Here q IS used as a shorthand notation for 
2-dimensional vector of the rescaled transverse momentum q = qt = pt/{l — z). The splitting 
function P{z,q,kt) is defined in eq.(|^) and the Sudakov form factor As^P^ax, zq) is given in 
eq.(|]). The first term in eq.(0) can be expressed as: 

A' = Go{x,k^)As{P^ax,ko) (18) 

with Go{x, k^) defined in eq.(H) and /cq being the starting scale for the evolution. Thus gives 
the probability of moving from xq, fcto and k^ to x, kt and Pmax without further branching. 

In a backward evolution, we start from gluon k^ (see Fig. [I|) and go successively down 
in the ladder until we end up at gluon k^. Therefore we need to reconstruct the branchings 
ki —>■ ki_i + qi, where fcj refers to the four vector of the t-channel gluon closest to the hard 
interaction, ki^i to the one closer to the proton and q^ to the gluon emitted between i and i — 
The momenta of such a branching with given Xj, ku and g^+i are reconstructed as follows:^ 



a. The value q^ = Ziqi is calculated from d?qi/qfAs{qi+i,qi) with the Sudakov form factor 
defined in eq.(|). Note, that due to the angular ordering constraint only the product z^gj 
can be obtained, whereas in the DGLAP case with transverse momentum ordering, 
could be obtained already here. 

h. Next the splitting variable zi = Xi/xi^i needs to be generated. The correct splitting 
function given in eq.(y) cannot be used, since it involves the variable q^ (in the argument 
of as), which is not known at this stage. Thus, Zi is generated first according to a 
approximate splitting function P^pp^\ 

P^ppr^z,, kl) = + ^il^ > (19) 

Zi i Zi 

which is always larger than the true splitting function. The limits on Zi are given by: 

^iQo 

Xi < Zi <1 . 

Also the z range needs to be larger than the true one: Zi < 1 — Qo/qi, but qi is not 
determined at this stage. Having Zi generated, the true qi can be calculated with qi = qi/ Zi. 
If now Zi > 1 — Qo/qi, the branching is rejected, and a new one is generated by going 
back to step a. 

■^I am very much grateful to G. Salam for suggesting this approach for a backward evolution. 



9 



c. Now the transverse momentum of gluon {i — 1) can be calculated, by generating a random 
angle around g^: kl_^ = fc| + g2 ^ 2y^^cos0. 

(i. A branching is accepted with the probability: 

^ 99 '^til f^max I J 

where Amax{xi) is the maximum value, A can take for any value x > Xi and any value of 

f^ti Qf 

This then completes the reconstruction of one branching. In the next branching generated the 
values of Xi-i, ku-i and qu play the role, the values Xi,kti and qu+i had before. The cascade 
ends, with the probability: 

A {xj-i, fcf^„i, qj ) ^^^-^ 
A{xi-i, qf) 

This means that a further branching would have come from the initial gluon distribution. 

Due to the complicated structure of the splitting function and A, this procedure is quite 
inefficient. It can be considerably improved by a different procedure in step c in the selection 
of the (f) angle and the transverse momentum kf^_i: the angle is now selected according to 
the probability: 

A qi ) 



A{xi-i, fcjmj^, qi ) 



(22) 



where fctmm — ^ti'^Qti^'^y^uyQu is the minimum transverse momentum allowed in this branch- 
ing. In step d the ratio of the parton densities (in eg. (^OD) is then changed correspondingly 

to: ^ ^ 

Xi~iA{xi-i, k^^j^^, q^ ) 

•^iAmax (-^i) 

This gives an improvement of a factor of ~ 2 in time. 

As mentioned already above, the true virtuality of the t-channel gluons can only be re- 
constructed after the full cascade has been generated. Now by going from the last (closest to 
the proton) gluon, which has virtuality k^ = kto/{l — xq), forward in the cascade to the hard 
scattering process, the true virtualities of kf are reconstructed. At the end, the gluon entering 
to the quark box will have a larger virtuality than without initial state cascade. Thus a check 
is performed, whether the production of the quarks is still kinematically allowed. If not, then 
the whole cascade is rejected, and the event without cascade is kept. 



4 Results 

In this section, results obtained from the hadron level Monte Carlo program Cascade, de- 
scribed in the previous chapter, are compared to those from Smmod. In all cases, no "consis- 
tency constraint" is apphed. The unintegrated gluon distribution was obtained from Smmod. 
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At parton level, both programs are expected to be identical, but differences can occur because 
of normalization problems, the finite grid size used to define of the unintegrated gluon density, 
and the effect of the virtuality of the gluon entering the hard scattering process, which is not 
properly known in the backward evolution approach. 




Figure 4: Comparison of the cross section obtained from the backward evolution Monte Carlo 
Cascade (solid line) with Smmod (dashed line) both at parton level only. The upper 3 plots 
show the cross section is shown as a function of the quark rapidity rjq, the quark transverse 
momentum pt and the transverse momentum of the quark pair p'^"'^^ . The lower 3 plots show 
the cross section as a function of the gluon transverse momentum kt and the multiplicity and 
transverse energy flow of the gluons from the initial state cascade are shown as a function of 
the rapidity rj 

In Fig. ^ the cross section as a function of the rapidity of the quarks rjg, the quark transverse 
momentum pj, the transverse momentum of the quark pair p^""^^ and the gluon transverse mo- 
mentum kt obtained from Cascade are compared to the ones obtained from Smmod (dashed 
line). A rather good agreement is obtained. Also shown is the multiplicity and the transverse 
energy flow as a function of rapidity for the gluons of the initial state cascade, and the pre- 
diction from Smmod for comparison. This shows that both the cross section and the initial 
state cascade are well reproduced within the backward evolution approach. However, the differ- 
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ences seen are due to uncertainties in the determination of the gluon density function. Further 
improvements needed here to optimize the grid and the extraction of the gluon density. 




In Fig. ^the structure function F2 obtained from Cascade is compared to the one obtained 
from Smmod (the same as the sohd hne in Fig. H). At large the results agree rather nicely, 
whereas at small differences are seen, due to uncertainties in the normalization of the 
gluon distribution. In Fig. |^ the prediction of the cross section for forward jet production is 
shown. The parton level results from CASCADE are similar to the ones obtained from Smmod, 
keeping in mind the differences already seen in F2. Also shown are the hadron level results 
from Cascade, which differ from the parton level results by ~ 20 %, which is known as the 
hadron level correction. 

In Fig. 1^ the transverse momentum distribution for single particles are shown in different 
bins of the rapidity. We observe a rather good description of the experimental measurements. 
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Figure 6: Comparison of the cross section for forward jet production obtained from the backward 
evolution Monte Carlo Cascade (dashed line) at parton level with Smmod (dotted line) at 
parton level. Also shown is the prediction from Cascade at hadron level (solid line), a. — c. 
The cross section for forward jet production as a function ofx, for different cuts inpt compared 
to HI data [^\] (a.~b.) and compared to ZEUS data [\^] (c). d. The cross section for forward 
jet production as a function of E^/Q"^ compared to f^^J. 



5 Conclusion 

In this article I have shown, that the CCFM evolution equation with a proper treatment of the 
non-Sudakov form factor, gives a good description of the structure function F2 and at the same 
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Figure 7: The transverse momentum distribution of particles in different bins of rapidity. The 
prediction of the Monte Carlo CASCADE at hadron level is compared to the measurement of 
HI 

time also describes the cross section for forward jet production. This resuhs were obtained 
from Smmod, a modified version of the forward evolution Monte Carlo program Smallx [ 
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I], 1^, where the important change was the modified treatment of the non-Sudakov form factor. 
Whereas for F2 equally good descriptions are obtained with and without the "consistency 
constraint" , the description of forward jet production improves if no "consistency constraint" 
is applied. 

To fully simulate the hadronic final state of small x processes including parton emissions 
according to the CCFM evolution equation, a new hadron level Monte Carlo event generator, 
Cascade, using a backward parton cascade approach, has been developed. The unintegrated 
gluon density as a function of x, k"^ and was obtained from Smmod. 

I have shown that the parton level results obtained from the forward evolution program 
Smmod agree rather well with those from the backward evolution program Cascade. The 
hadron level results obtained from Cascade show rather nice agreement with measurements of 
the HERA experiments, that are sensitive to small x parton dynamics. With standard hadron 
level Monte Carlo programs using DGLAP type parton evolution and direct photon interactions 
only, a similar level of agreement could not be achieved. 
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